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Abstract 

We extend the formulation of relativistic Coulomb sum rules to account for the 
average effects of nuclear binding on the initial and final states of ejected nucleons. 
Relativistic interactions are included by using a Dirac representation adapted from a 
vector-scalar field theory. The scalar field reduces the effective nucleon mass M* and 
increases the relativistic effects of recoil and Fermi motion. We consider two models 
for the off-shell behavior of the nuclear electromagnetic current, and demonstrate that 
the sum rule is accurate for applications to data over the interesting range of M* and 
three momentum q. We further indicate that the form of the sum rule is sufficiently 
general to accommodate a broad class of off-shell form factor models. 
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1 Introduction 



The subject of this paper is the Coulomb sum rule (CSR) for inelastic electron scattering [Q, 
in a reformulation which extends its validity to relativistic momentum transfers q M, 
where q = |q| is the three-momentum transfer to the target nucleus, and M is the nucleon 
mass.Q The purpose of this work is to provide practical methods for the analysis of the 
longitudinal (e, e') response of nuclear targets, at the higher energies now available at the 
Continuous Electron Beam Accelerator Facility (CEBAF) and other electron accelerators. In 
a previous paper 0, we derived a relativistic Coulomb sum rule (RCSR) which incorporates 
the relativistic effects of nucleon recoil at large g, as well as of Fermi motion. In the present 
work we extend this theoretical approach to include the effects of nuclear interactions on the 
RCSR. 

In our recent article, we discussed in detail the assumptions under which a RCSR could 
be derived. The basic approximation is that only nucleon, i.e., as opposed to antinucleon, 
degrees of freedom enter the Coulomb response in the spacelike regime accessible by (e, e') 
experiments. We call this the "nucleons-only" approximation; it ignores effects of antinu- 
cleons, but includes fully the relativity of the nucleons. Following conventional treatments, 
we adopt an impulse approximation which ignores explicit contributions from the exchange 
of charged mesons. We include nucleon anomalous moments and elastic form factors, which 
are important at higher energies. The final assumption is less conventional, but absolutely 
necessary. In order to derive a non- energy-weighted sum. rule, it must be possible to factor all 
dependence on the photon energy uo from the current matrix element, otherwise dynamical 
effects enter into the Coulomb response function and complicate the isolation of correlation 
effects. In this paper, we restrict our attention to form factor models which have factorable 
dependence on u. 

We define the nuclear Coulomb sum in terms of the Coulomb response function Wc'{uj, q) 
and the proton electric form factor Ge,p{Q'^)'- 

where the lower limit cu^ excludes the quasielastic peak, and the upper limit q restricts the 
integration to spacelike four-momenta. Under the assumptions of Ref. mentioned above, 
we obtain a RCSR which can be expressed in terms of one- and two-body contributions: 

S(q)^S«(q) + C(q) + S(2„)(q), (1.2) 



^For a derivation of the nonrelativistic CSR in second quantization, see Ref. 
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where the one-body contribution is of the form 



EW(q)^25:n<,(p)r.(p,q), (1.3) 

pa 

C(q) includes two-body correlation information in momentum space, and S|j^)(q) is the 
uncorrelated two-body part which is related to the square of the nuclear elastic form factor, 
as discussed in Ref. [|]. 

In the one-body term (|1.3|) , ^^(p) is the nucleon momentum distribution function for 
isospin projection a and one spin projection, and ro-(p, q) is a kinematic factor which arises 
due to relativistic nucleon recoil and Fermi motion in the target. In the nonrelativistic limit 
(g « M), we have Vp—^l and r„ 0, which leads to S^^)(q) Z; then the sum rule ( |1.2| ) 
is the result of Ref. The relativistic effects are all in the functions ro-(p, q), representing 
recoil of the struck nucleon and Fermi motion in the target ground state. In Ref. [0] we 
further showed that ( |1.3| ) does not depend strongly on the details of ^^(p), but only on 
the lowest momentum moments, e.g., (p^)(j- This leads to a method of evaluating (|1.3| ) 
accurately in a weakly model-dependent manner, and in principle permits the extraction of 
the correlation function C(q) from the the experimentally measured Coulomb sum ( |1.1|) , 
using (|r|). 

In the present paper, we investigate what changes are required in the sum rule when 
the Coulomb sum S(q) is modified by relativistic two-body interactions. We study the 
effect of the average interaction in the nucleus using the mean field approach of quantum 
hadrodynamics, a relativistic field theory for nuclear physics 0]. In particular, we consider 
QHD-I, which includes vector and scalar isoscalar mesons only. We then reformulate the 



Coulomb sum rule of Ref. [pi so that results (1.1)-(|1.3|) have a similar structure, but with 



modifications reflecting the mean-field effects of these relativistic interactions. 

These modifications have two main effects on the RCSR: First, there are kinematical 
effects resulting from mean-field interactions of the nucleons, which are represented in QHD- 
I by a reduction in the effective nucleon mass M* in the medium. Consequences of a reduced 
effective nucleon mass M* for the Coulomb sum have been considered previously in a Fermi 
gas model Q . These calculations include interactions of both initial and final plane- wave 
nucleon states with the mean fields, through the effective mass M*, which can be interpreted 
as binding in the initial state, and final state interactions of the ejected nucleon with the 
nucleus. Chinn, Picklesimer and Van Orden ||^, || have studied the effects of final state 
interactions on the Coulomb response of a Fermi gas, using more realistic interactions, and 
have seen similar effects to those seen due to M* /M < 1. Second, the electromagnetic 
coupling of the nucleon in medium may be modified by the mean fields, entering through the 



3 



off-shell behavior of the nucleoli elastic form factors. Both modifications introduce a degree of 
model-dependence in the RCSR which is not present in the nonrelativistic formulation, nor in 
the relativistic formulation of Ref. 0. We show how these features can be incorporated into 
the theory to allow the evaluation of the one-body contribution S^^''(q), and the subsequent 
extraction of the two-body correlation function C(q) from the measured Coulomb response. 

This paper is organized as follows: In Section 2 we introduce the basic formalism to 
include relativistic mean-field effects in the Coulomb sum rule (RCSR). We introduce two 
models (F and G) for the electromagnetic charge operator, and investigate the resulting 
behavior connected to different off-shell assumptions for the nucleon elastic form factors. In 
Section 3, we derive a modified version of the RCSR ( |l.l|) -( pr3D , concentrating on the explicit 
changes to the one-body term S'^^^(q) in off-shell Models F and G. In Section 4, we illustrate 
the operation of the sum rule in a simple nuclear system: uniform nuclear matter treated in 
the mean-field approximation, with nuclear binding effects incorporated using QHD-I. We 
examine the sensitivity of the RCSR to M* and to the choice of off-shell models (F and G), 
focussing on the convergence of the moment expansion in each case. We further demonstrate 
that the particular form of the RCSR given here is applicable to both models, and argue that 
the same form should also be valid for a broad class of form factor models. In Section 5, 
we draw conclusions, give guidelines for the application of the RCSR to data, and indicate 
important directions for future work. 



2 Formalism 

In this section, we review the formalism for electron scattering from nuclei, as it pertains 
to our development of relativistic Coulomb sum rules. We first give some standard results 
of the plane- wave impulse approximation (PWIA), which is used for the analysis of (e, e') 
experiments on nuclei, in a single-particle basis which accounts for the vector and scalar 
interactions of QHD-I. We then give two models for the nucleon form factors off shell, which 
will be used in the next section to illustrate the sensitivity of the sum rule to different off-shell 
assumptions. 



2.1 PWIA in M* basis 

We begin with the differential cross section for the scattering of ultrarelativistic electrons 
from nuclear targets, which is commonly written in the form 



da. 



dVt'dE' dVt' 



q4 



i^+tan^^)iyH^,q) 



(2.1) 
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where = (cu, q) is the four- momentum transferred from the electron to the nucleus via 
virtual photon exchange, and = |qp — cj^ > 0. The longitudinal contribution in ( |2.1| ) has 
been expressed in terms of the Coulomb response function 

W^c(^,q) = T.MUqM? -Ef + E,), (2.2) 
/ 

where \i) and |/) denote initial and final nuclear many-body states, respectively. 

In general, the electromagnetic current density operator Jfj,{q) may include contributions 
from both nucleons and charged mesons. Meson exchange current (MEC) contributions have 
been considered by Schiavilla et. al. Q, for example, but are not included here. Including 
only nucleons with electromagnetic form factors, the current density operator can be written 
in the form 

J^(g) = J e'^-^ ^(x)r^(g)z^(x), (2.3) 

where V'(x) is the (Schrodinger picture) field operator for a point Dirac nucleon, and 
represents the electromagnetic coupling at the 'yNN vertex. There is also a sum over proton 
and neutron isospin projections, which will be suppressed until needed. Any effects of charged 
mesons not included in the current operator ( ^.31 ) may be interpreted as two-body effects, 
as discussed in Ref. [0. 

To proceed to a relativistic sum rule, we shall expand the field operators ^(x) of (|2.3| ) 
in a plane wave basis, i.e., in a PWIA. In Ref. [Q], we took this to be the free basis, i.e., 
the momentum eigenstates of the free Dirac equation. In the present paper, we shall instead 
expand the field operators in a basis of plane waves moving in the presence of uniform 
(isoscalar) Lorentz scalar and vector potentials, such as would be generated by scalar and 
vector mesons in a relativistic field theory of nuclei, e.g., QHD-I. With this modification, 
the theoretical development is formally similar to that of Ref. [0]. The consequences of 
the change of basis show up after the nucleons-only approximation, in which antinucleon 
degrees of freedom are removed approximately from the spacelike Coulomb response. In the 
following, nucleons-only will refer to nucleons in the presence of scalar and vector potentials 
in the target. The effect of this modification is to include efficiently the effects of these 
mean-field potentials in both the initial and final target states. 

We therefore use the plane-wave solutions of the Dirac equation 

r{id^ - g.V,) - (M - gs^)U{^) = 0, (2.4) 



^For notational simplicity, we assume a nondegenerate target ground state; the results are easily gener- 
alized to unpolarized targets with J^O. 
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where the scalar and vector potentials are written in the forms gs(f) and (^t^V^, respectively, 
corresponding to the scalar and vector fields (j) and V^, and the associated coupling constants. 
Here, as in the mean-field solution for uniform nuclear matter at rest, V^ = 5|JVo, and Vq and 
are constants. The solutions to ( p.4|) are discussed in detail in Ref. In the case of 
uniform fields, the energy eigenvalues of ( p.4| ) take the simple form = g^Vo ± -E*, where 
= \J p^ + M*^ and the effective nucleon mass is defined M* = M—gs4>. With the nucleons- 
only approximation, we need only consider positive-energy solutions to ( p.4|) . These are 
plane wave solutions of momentum p, which obey the equation 



7°S;-7-p-M*]w,(p) = 0. (2.5) 

Explicit forms for the interacting solutions Us{p) can be obtained directly from the free 
solutions, given in Appendix A of Ref. 0, by making the replacement M— i>M*. Since the 
vector potential g^Vo appears additively in the nucleon eigenenergy, it does not appear in 
( p.5|) or its plane-wave spinor solutions. 



In Ref. Q we argued that the Coulomb response function (|2.2| ) for spacelike (cj < |q|) pho- 
ton exchange is dominated by nucleon (NN) contributions to the current matrix elements, 
and that antinucleon {NN) and pair {NN) terms could be neglected. This is an exact result 
for a uniform free Fermi gas, and leads to the nucleons-only approximation for interacting 
nuclear systems. The presence of a strong scalar field in the nucleus induces mixing of free 
and A" states, particularly for nucleons of high momentum, as in the final states of ( |2.2| ). 
However, transforming to the plane- wave basis formed from the solutions of ( |2.5| ), hereafter 
referred to as the "M*-basis" , removes this mixing by the potentials. The nucleons-only ap- 
proximation is again adopted for interacting nuclei, but here refers to nucleons of mass M*. 
This use of the M*-basis for final states implies substantial interaction of excited (ejected) 
nucleons before leaving the target, and is probably a better assumption for large nuclei than 
for small. Then the electromagnetic current operator takes the formQ 

(2.6) 

P ss' y^-C/p+q Y P 

where a|,g and is a creation operator for a nucleon with spin projection s and momentum p. 



The energy denominators in (|2.6|) refiect the normalization of the plane-wave spinors to 2Ep 
particles/volume. The formal derivation of the RCSR now follows closely that of Ref. [Q, 
once we have discussed the form of the electromagnetic vertex operator T^{q). 

•^For simplicity, we discretize the sum over the momentum p. 
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2.2 OfF-shell nucleon form factors 

It is conventional to express the 'jNN vertex operator in the form 

r,{q) = F^^^ + z^F,a,,q\ (2.7) 

where k is the nucleon anomalous magnetic moment, M is the free nucleon mass, and Fi 
and F2 are the Dirac and anomalous form factors, respectively. In general, Fi and F2 are 
scalar functions of p, p' and q. This form is sufficiently general for matrix elements between 
nucleon states (antinucleons excluded) as in the nucleons-only PWIA. Since the tensor a^u 
is antisymmetric, only the three- momentum q enters explicitly the Coulomb operator Tq. 
Any remaining dependence the photon energy u enters Fq only through the form factors Fi 
and F2. 

For scattering from a free nucleon, the form factors Fi and F2 depend only on the 
scalar Q^. To ensure the correct charge and magnetic moments for free nucleons, they 
are normalized at = according to -Fip(O) = F2p{0) = -F2„(0) = 1 and -Fi„(0) = 0. The 
form factors Fi and F2 are obtained from (e, e') scattering data, usually in terms of the more 
convenient Sachs electric and magnetic form factors Ge and Gm'- 



«i^2(Q') 



GEiQ') + TGMiQ') 
1 + r 

(2.8) 

Gm{Q')-Ge{Q'') 
1 + r 



where r = Q /AM , and Q = —{k^ — k'^) , i.e., the photon four-momentum is determined 
by the four-momentum transfer at the jee vertex. For a single nucleon in free space, both 
energy and momentum are conserved at the 'yNN vertex as well, and we have p'f^=Pij, + q^- 
See Appendix A for an alternative form of (2.7) which is expressed in terms of the Sachs 
form factors. We use the standard parameterization of the Sachs form factors, along with the 
assumption G_Bn(Q^) = 0, as discussed in Ref. [Q. This choice has the convenient feature that 
all Sachs form factors are proportional, which satisfies a condition assumed in our derivation 
of a non- energy-weighted sum rule in Ref. 

For interacting nucleons in a nucleus, one needs information about the form factors 
off their free mass shell, which is simply not known. This leaves considerable freedom to 
extrapolate off shell from the known on-shell forms. Given no other information, a common 
approach has been to assume that the dependence on of the form factors Fi((5^) and 
F2{Q^) does not change off shell from that for a free nucleon. This implies, at least in part, a 
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separation of the dynamics of the electromagnetic structure of the nucleon from that of the 
nucleus. However, even this assumption does not uniquely specify the off-shell form factors, 
because of the freedom to transform r^(g) using the Gordon identity to other operators with 
other form factors, each of which are equivalent on shell, but not off shell. This issue has been 
considered by de Forest |jlO|], and by Chinn and Picklesimer who have given examples 
of possible choices of extrapolation, and studied the sensitivity of the response functions 
and of the Coulomb sum to those choices. Our method is similar in principle, but is based 
specifically on the Gordon transformation which relates Fi and F2 to the free Sachs form 
factors Ge and Gm^ given in ( p.8|) for free nucleons. We consider two possibilities: that the 
functional dependence on Q"^ of the functions F2{Q'^) is unchanged off shell (Model 

F), or alternatively, that the functional dependence on of the functions Ge{Q'^), Gm{Q^) 
is unchanged off shell (Model G). Other choices are possible. For example, in principle, the 
vector field Vq may also enter the current operator. We restrict our attention to Models F 
and G, which illustrate the important issues and do not depend on Vq. 

The first choice. Model F, is the most common; here it is simply assumed that the Dirac 
and anomalous form factors Fi{Q'^) and F2{Q'^) are unchanged in the nuclear environment, 
i.e., that they are given by ( |2.8| ) where both the Sachs form factors and the kinematic 
variable r = Q'^/AM'^ are evaluated at the actual momentum transfer of the experiment, 
i.e., at the momentum = k^ — k'^ found at the electron vertex. This is in some sense 
a minimal assumption, in that the effects of the scalar field enter explicitly the current 
operator J^{q) only through the nucleon field operators. In spite of its apparent simplicity, 
this model has the peculiar feature that the Dirac and anomalous magnetic moments are 
treated differently in the nuclear medium. This follows since the Dirac moment (e/2M*) 



scales with the nucleon mass M*, as can be seen from the solution of ( p.4|) for a nucleon at 
rest in the presence of a uniform magnetic field, while the anomalous moment {Ke/2M) is 
unchanged in the medium. This behavior is possible for a theory with point Dirac nucleons 
dressed by charged meson fields, such as QHD-11 since the anomalous moments in such 
a theory are typically included explicitly by hand, following a solution of the many-body 
problem using point nucleons. However, this behavior is not likely for a theory with internal 
nucleon electromagnetic structure, such as QCD, since there one expects that the Dirac and 
anomalous magnetic moments may have similar origins, and therefore may respond in a 
similar manner to the nuclear environment. 

An alternative assumption, which we call Model G, is that the Sachs electric and magnetic 
form factors Ge{Q'^) and Gm{Q'^) are unchanged in the nuclear medium. To express r^(g) 
in the form ( |2.7| ), the transformation ( |2.8| ) is now performed in the M*-basis: the off-shell 
forms of Fi and F2 are now given in terms of the free Sachs form factors by 
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Ge{Q') + t*Gm{Q^) 
1 + f* 

(2.9) 

M Gm{Q')~Ge{Q') 



M* 



where the new variable f* = [q2-(E;^q-E;)2]/4M*2. In contrast to Model F described above, 
the off-shell choice ( |2.9| ) leads to a total nucleon magnetic moment equal to + e/2M*. 
In particular, the effective Dirac and anomalous magnetic moments behave similarly in the 
nuclear medium. (Note that the factor M/M*, appearing on the right-hand side of F2 in 
( p.9| ), has the effect of replacing k/2M by k,/2M* in ( p.7| ).) Also, since f* depends only on 
the momenta p and q, the photon energy uj now enters the form factors Fi and F2 only 
through the Sachs form factors Ge^Q"^) and Gm^Q"^). This can be seen by expressing Fq in 
the alternate form given in (A5) of Appendix A. In the next section, we will see that Model 
G leads to relativistic Coulomb sum rule which is the direct analog of that derived in the 
free PWIA, but with M^M*. 



3 Coulomb Sum Rules 

The derivation of a RCSR based on the interacting PWIA of Section 2 is formally similar to 
that based on the free PWIA of Ref. We begin with the Coulomb response function ( |2.2| ), 
and formally evaluate the spacelike Coulomb sum ( |1 . 1| ) by integration over the photon energy 
oj. As described in Ref. to arrive at a non-energy-weighted sum rule it must be possible to 
factor all dependence on u (which here enters through Q"^) from the current matrix element 
in (|2.2|). This can be accomplished by requiring that the ratio Fo(g)/G'£;p(<5^) be independent 
of UJ, since the plane-wave spinors appearing in (|2.3| ) are functions only of the 3-momenta p 
and q. For a system of Dirac protons, for which Fo = 7o, this is satisfied trivially. In Model 
G, it is satisfied by explicit construction, using the assumption of proportional Sachs form 
factors. For more complicated off-shell models, such as Model F, which do involve explicit 
dependence on u, the derivation of a RCSR may require further assumptions. In this article, 
we will assume that the above condition is satisfied, and outline the remaining steps which 
lead to a RCSR. 

We next use closure to perform the sum over final states in the squared matrix element. 
In the nonrelativistic Coulomb sum rule, in which the integration is over all u, the use of 
closure here is exact. In the relativistic case, the use of closure over the spacelike states alone 
requires certain assumptions about the spacelike nuclear excitation spectrum, as discussed 
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in Ref. 0]. These arguments were based primarily on the example of a Fermi gas. We 
make similar arguments here, and assume that the spacelike spectrum is saturated by the 
nucleons-only response, where now "nucleons-only" refers to positive-energy baryons of mass 
M* . After performing the sum over final states |/), the momentum-space anticommutation 



relations are used to separate one- and two-body terms, as in (|1.2| ). The relativistic recoil 
function which appears in (|1.3|) is of the form 



2 



'^^(p, q) = ^ hs'sA'P^ q) ^ (3-1) 



2 

where the matrix element is defined by 

and 'Uso-(p) represents a nucleon with momentum p, spin s and isospin projection a. The 
above expressions are identical in form to those obtained in the free PWIA of Ref. 0, but 
with the replacement M^M*. 

Although the one-body term, evaluated using ( |3.1| ) and ( |3.2| ) would account exactly for the 
relativistic effects of nucleon recoil and Fermi motion, it requires knowledge of the nucleon 
momentum distribution ?2o-(p) in the target ground state, which is not normally known 
precisely. We therefore expand the recoil factor r^ip, q) in powers of the nucleon momentum 
p, as in Ref. ||^. Assuming spherical symmetry, the one-body term takes the forniQ 

S«(q)=E N^nM) (p'). , (3.3) 
where the momentum moments are defined 

(p^)<^ = i|- E ^'^(p) p^ (3-4) 

and Nf^ = Z,N for a = p,n, respectively. In Ref. 0, we explained that only the first few 
moments are required to obtain an accurate result for the one-body term. In the next 
section, we will study how the convergence properties of this moment expansion depend on 
the parameter M*, and on the particular choice of off-shell model. We now give the precise 
forms of the recoil factor ricr{q) in Models F and G. 

''We have changed notation shghtly from Ref. in naming the expansion coefficients ria{q)- 
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3.1 Sum Rules in Model G 



We first consider Model G, in which we assume that the medium dependence of the nucleon 
electromagnetic coupling is described by the off-shell form factors ( ^.91 ) . By construction the 
form factors Fi and F2 depend on only through their linear dependence on the free Sachs 
form factors, which are assumed to be proportional. This ensures that the ratio '^q/Gep{Q'^) 
is independent of the photon energy uj. With this condition satisfied, and since the baryon 
spinors are also independent of cu, the matrix element js's,cr(p, q) is a function only of three- 
momenta and the separation of one- and two-body terms in S(q) proceeds as described 
above, leading to ( |1.3|) . 



Using the form factors ( |2.9| ) in ( |2.71) to evaluate the matrix element js's,(j(P; ^) we obtain 



q) = 4 ^Eiv, q) + '^m(p, q), (3.5) 
where we have defined the electric (E) and magnetic (M) recoil functions: 

Mp,q) = ^^^^^, (3.6) 

and 



'^M(p,q) - 



1 + f * 4E;£;;+q 



|p X q 



(3.7) 

2 



The functional forms of te and tm can alternatively be obtained by writing of (2.7) in 
terms of G^; and Gm, as explained in Appendix A. In ( |3.5|) we have also defined the nucleon 
electric charge e„ (in units of the proton charge) and the total magnetic moment fi^j (in 
nuclear magnetons): 

GEaiQ^) _ [ 1, foTa=p 
" GEpiQ^) lO, fora = n 

(3.8) 



GmAQ'^) _ ( 1 + K,p, for a = p 
Gep{Q^) I f^n, for a = n. 



The functional form of the recoil function (|3.5D is identical to that obtained in the free PWIA 



of Ref. 0, but with M— i>M*. This was ensured by choosing to include the factor M/M* in 
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( p.9| ). As described above, we expand ro-(p, q) in powers of the nucleon momentum p, which 
with (|1.3|), leads to an expression for the one-body term S(^)(q), as in (p.3[). The expansion 
coefficients are given in Appendix B, and are identical to those obtained in the free PWIA of 
Ref. 0], but with M— >M*. In this model, therefore, we can include the relativistic effects of 
a strong scalar field, in a straightforward extension of the RCSR derived in the free PWIA. 



3.2 Sum Rules in Model F 

We now return to Model F, in which we assume that the nucleon form factors are given 
by (|2.8| ). Now the photon energy uj enters not only through Q^, which appears in the free 
Sachs form factors, but also through r = Q^/4M^, which enters in the transformation (|2.8|) . 
This additional dependence on uj can not be removed by simply dividing by an overall factor 
Gep{Q^), as in Model G, without some further approximation. 

In order to proceed, we will first calculate the matrix element appearing in (|3.1| ), using the 
vertex operator 



and the form factors (2^). Since division by GspiQ'^) is not sufficient 
to render the matrix element a function only of the momenta p and q, this function will not 
yet lead to the one-body term of the Coulomb sum rule. However, it does serve as a useful 
starting point to illustrate the important issues. We obtain 



(3.9) 



where we have emphasized explicitly the dependence on Q^. The form of ( p.9|) is similar to 
(3.5) for Model G, but with "effective" Sachs electric and magnetic form factors, defined by: 



M* 



(3.10) 



GmAQ 



M* 



which reduce to the free Sachs form factors if M* = M. Inserting FilQ"^), F2{Q'^) of ( 12!S 
into (|3.1CI| ), we have 



G*EAQ',r* ^ 
Gep{Q^) 



GlsAQ',r*] 
Gep{Q^) 



1 I Mlf*n 
1 + T 



1 - 



M* -\ 

M 



l + r 



(3.11) 



r + 



M* - 

M 



l + r 
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which illustrates of "mixing" of (free) electric and magnetic contributions in this model. 
This mixing enters as a result of the reduced effective nucleon mass, and is distinct from the 
usual Lorentz mixing which occurs for a moving particle. Effective form factors of the sort 
(3.11) can also be seen in Ref. for example. 

To continue the derivation of a non- energy-weighted sum rule, we must make some as- 
sumption about the remaining ci;-dependence, which arises through the variable r in (3.11). 
A simple method is based on the excitation energy of a uniform Fermi gas in QHD-I, which 
is expressible simply in terms of the three-momenta p and q: u; = i?p_,_q — i?* . We therefore 
make the replacement in (3.11): 



to obtain energy-independent effective form factors. The factor (M*/M)^ arises directly from 
the definition of r following (2.8). With the replacement ( |3.12| ) in (3.11), the recoil function 
(3.9) is now independent of the photon energy tu, and the derivation of a non-energy-weighted 
Coulomb sum rule for Model F proceeds as for Model G. 

At this point, it is possible to expand (3.9) about p = and obtain a moment expansion 
for S*^^)(q), as in ( p.3|) . However, the resulting coefficients rjo-(q) are much more complicated 
than those which arise from (3.5) in Model G, due to the dependence on p of the effective 
form factors (3.11). In addition, we will see that the essential physics of Fermi motion enters 
( p.9|) only through the electric and magnetic recoil functions r£;(p,q) and ?"m(p, q), and 
can be ignored in the effective form factors themselves. It is therefore useful to consider 
an approximation scheme which will allow a more efficient evaluation of the one-body term: 
setting p = in (3.11) after applying (3.12), we can write 




(3.12) 



r^(p, q) = e;^(q) rs(p, q) + /^^^(q) rM(p, q), 



(3.13) 



where we have defined the "effective" nucleon charges and magnetic moments: 




p=0 



2M2 + M(E* - M*) 



{M* - M){E^- M*) 



2M2 + M*(E*-M*) 



2M2 + M*(E* -M*) 



(3.14) 




p=0 



13 



2M(M - M* 



2MM* + M*{E^- M* 
2M2 + M*(E*-M*) 



2M2 + M*(E*-M* 

following ( p.8| ). We will refer to expression ( p.l3| ) as the "factored moment expansion." The 
effective form factors now enter ( p.l3| ) as functions only of the momentum transfer q, but 
do not complicate its expansion in p. That this simplification is a valid approximation will 
be shown numerically in the next section. 



4 Sum Rules for Relativistic Nuclear Matter 

We now apply the sum rules derived in the previous section to a model nuclear system: 
uniform nuclear matter in the relativistic mean field approximation, as given by QHD-1. 
For this system the effects of nuclear structure on the Coulomb sum S(q) enter through 
the effective nucleon mass M* and through the nuclear electromagnetic current, as we have 
discussed in Section 3. We are interested in the dependence of the Coulomb sum on M* and 
on the choice of off-shell model (F or G) for the current. We then demonstrate how to apply 
the sum rule methods of Section 3 to the model system, as if it were measured Coulomb 
response data. The analysis will necessarily be model dependent, through the choice of 
M* and the off-shell model (F or G), as well as through the moments (p*)o- of the nucleon 
momenta in the target. We shall examine the accuracy of the moment expansion for the 
"best-case" analysis, that is, for which the sum rule parameters match those of the assumed 
nuclear model. This will illustrate the convergence properties of the moment expansion in 
each model, and the validity of the factored moment expansion in Model F. 



4.1 Test in Fermi Gas Model 

We begin by computing the Coulomb sum S(q) for a uniform Fermi gas of nucleons moving 
in uniform vector and scalar potentials. This has been studied previously |^ . The result 
can be written in the form 

7^ II^a(p)[l-?^a(p + q)] |js's,<x(P,q) , (4.1) 

where the matrix element js's,a{p, q) is defined in ( |3.2| ). We will evaluate ( |4.1| ) by numerical 
integration over p and 6, the angle between p and q. The result is independent of Vq, for 
reasons given following ( |2.5| ). The factor ?T,o-(p) =0{pFcr — \p\) restricts the sum over initial 
states to include only those states which are occupied in the nuclear ground state, and the 
factor [1— ?T,o-(p+q)] ensures that final states for which |p+q| <PFa are excluded. This Pauli 
exclusion is the only source of two-particle correlations in this simple model. 
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We shall consider two particular values of M* which arise in QHD-I, as well as M* = M. 
In the mean-field theory (MFT), static polarization of the baryon vacuum in the nuclear 
ground state (due to the scalar field 0) is ignored in the coupled field equations. In the 
relativistic Hartree approximation (RHA), static vacuum polarization is included at the one- 
loop level.0 The particular values of M* are obtained from the self-consistent solutions 
for uniform nuclear matter at saturation. Taking Z = N and using the input parameters 
Pf = 1.42 fm"^ and BE/N = -15.75 MeV/nucleon, leads to the values M*/M = 0.556 in 
MFT, and M*/M = 0.718 in RHA. Although it may be argued that the model assumptions 
of MFT are technically more consistent with the nucleons-only approximation, it may also 
be true that the resulting effective mass M* is too small due to the simplicity of the model, 
as compared to optical potential phenomenology (see, e.g., Ref. fl^). In the calculations to 



follow, we will show results using both the MFT and RHA values of M* to illustrate the 
kinds of effects which can be expected in this range. 

For applications to data, it is convenient to cast the RCSR in a slightly different form, so 
that the saturation of the sum rule is more apparent. As in Ref. [Q, we define the "modified" 
Coulomb sum Q 

S4q)^^^ A = I, II, II I... (4.2) 



where the recoil "correction" factors are defined 



N 



rop + —ron 



riiiq) = ri(q) + 



(4.3) 



(4.4) 



and so on for higher orders. [] The expansion coefficients rjo-(q) are given in Appendix B. 
The moments (p*)fT are defined in (3.4). For applications of the RCSR to S(q) of ( [4.1|) , we 
will assume the sharp distribution ncr{p) = 0{pFcr — |p|) for the evaluation of the momentum 
moments. Integrating this distribution over the Fermi sphere of radius ppa^ leads to the i^^ 
moment, given by 

(P% = ^ • (4.5) 



This is consistent with the PWIA and the evaluation of the RCSR on ( [4.1|) . For applications 
to actual data, however, it would be preferable to use the most accurate values available for 
these moments, e.g., using experimentally determined ?^cr(p) from {e,e'p). 

5 "MFT" and "RHA" are used here in the conventional sense defined in Ref. ji) . 
^Note the factor 1/Z, which is not present in Sa as defined in Ref. 
^In this notational scheme, rA(q) includes terms through 0{p'^^^^). 
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Following ( p..2| ), the above definitions lead to an approximate sum rule of the form 



^^(q) ^ 1 + 



1 



C(q)+S(i)(q) 



(4.6) 



as discussed in Section 5 of Ref . . The Fermi gas model includes only Pauli correlations, 
which vanish identically for q>2pp. Consequently, an application of the RCSR to ( [4.1| ) will 
be considered "accurate" if the modified Coulomb sum S'^(q) tends to unity for momenta 
beyond the range of these correlations, i.e., if 5'yi(q) — > 1 for q > 2pp. This is frequently 
referred to in the literature as "saturation" of the sum rule. 



4.2 Numerical Results 

In the following examples, we take Z = N. To demonstrate the importance of the effects of a 
reduced effective mass M*, we first consider a system of Dirac protons, for which ep = fip = l 
and en = fJ'n = 0. In Ref. we showed for M*/M = 1 that the lowest-order sum rule, 
Si{q), is accurate to within ~ 1% for this case (see Figure 2 in Ref. 0). In Figure 1, the 
dotted and dot-dashed curves show the Coulomb sum S(q) calculated with M*/M = l and 
M*/M = 0.556, respectively. In this example, we see a 10-20% reduction in the Coulomb 
sum over this momentum range, due to the reduced effective mass M*. This effect has been 
noted previously for lower momenta P, |^. The dashed and solid curves show the lowest 
order sum rule, 5'/(q), applied to S(q, M*) in two different ways: In the dashed curve, the 
recoil factor r/(q) has been evaluated at the free nucleon mass M, while the solid curve shows 
a consistent application of the sum rule, in which r/(q) has been evaluated at the reduced 
mass M*. The dashed curve demonstrates that if the effective nucleon mass is substantially 
reduced in the medium, then one must account for this in the application of the RCSR. 
The solid curve demonstrates that even when the effective mass is greatly reduced, e.g., 
M* /M = 0.556, the lowest-order sum rule 5'/(q) applied to a system of Dirac nucleons is 
accurate to within ~1%, if the appropriate value of M* is used. 

In Figure 2, we show the RCSR applied to a uniform system of nucleons with anomalous 
magnetic moments, for M*/M = 1. (Here Models F and G are identical.) A similar figure 
was shown in Ref. 0, but it is useful to review the main results here for comparison. The 
unmodified Coulomb sum S(q) is shown by the dotted curve, and indicates an enhancement 
relative to the free Dirac result. This is a consequence of the nucleon anomalous moments, 
as explained in Ref. p. The lowest-order sum rule, Si{q), overshoots the "correct" result, 
i.e., saturation, increasingly with higher momentum, and indicates the importance of Fermi 
motion effects when anomalous moments are included. In Ref. we explained that this is 
because the magnetic contribution is enhanced relative to the electric by a factor fip+fJ.^ — 11-4. 
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This prevents the cancellation of electric and magnetic Fermi motion effects, which occurs 
nearly exactly for Dirac nucleons. The striking feature is that the O(p^) result, Sii{q) is 
accurate to within ~ 1% over this momentum range. We also show the O(p^) result, ^'///(q), 
which is accurate to better than 0.1% over the same range. This we interpret as numerical 
evidence of convergence of the moment expansion in this model. The results shown here will 
be modified for M*/M < 1, and will also depend on the choice of form factor model (F or G). 
We are interested in the accuracy of sum rules designed to account for these modifications. 

We now consider M*/M < 1. In Figures 3a and 3b we show results for the RCSR in 
Model G for RHA and MFT values of M*, respectively. As for Figure 2, the Coulomb E(q) 
is enhanced relative to that for Dirac nucleons, due to Fermi motion effects, when anomalous 
moments are included. This model also has the rather pecuhar feature that for M*/M~ 1/2, 
we have E(q)~l for q>2pF. This value is coincidental (see Sec. 4.3 for further discussion). 
For still lower values of M* the Coulomb sum can be significantly greater than unity. This 
is contrary to the common notion that relativistic effects necessarily suppress the Coulomb 
sum relative to unity. Note that the moment expansion converges more slowly for smaller 
M*, as can be seen by comparing Sii{q) in Figures 3a and 3b. Still, the RCSR to O(p^) is 
reasonably accurate, even for M*/M~ 1/2. The very high accuracy of Siii{q) is interpreted 
as numerical evidence of convergence in this model. Indeed, for q>2pF, Sm — l.OO for the 
RHA mass and Sm ^1.01 for the MFT mass, at g = 4 GeV. 

In Figures 4a and 4b we show analogous results for a nuclear system under the off- 
shell assumption of Model F. The Coulomb sum E(q) behaves qualitatively differently in 
this model, with a minimum developing around g ~ 2 GeV as M* decreases to the MFT 
value. The differences originate from the additional ^-dependence in the effective form 
factors G*f^/Gep and G*j^^/Gep, which are constants in Model G. As in Model G, Si{q) 
is not adequate, and illustrates the importance of Fermi motion effects when anomalous 
moments are included. Also as for Model G, 5'//(q), shown by the dashed curve, is reasonably 
accurate over this range of momenta and M*. The improvement in Siii{q) is interpreted as 
numerical evidence of convergence of the moment expansion. The most important feature 
of the sum rule results just presented for Model F is the success of the factored moment 
expansion, which allows the use of the same expansion coefficients as in Model G, i.e., those 
given in Appendix B. This approximation neglects only the Fermi motion effects which 
enter through the ratios [G*/GepY in (3.9). To illustrate its accuracy, we have evaluated 
Sji{q) in a "consistent" moment expansion which correctly includes to C(p^) additional 
Fermi motion effects entering through these ratios. This result is shown by a dot-dashed 
curve, which is nearly identical to the "factored" result, shown by the dashed curve. These 
are small corrections on the scale of the other kinematical effects, lending support to the 
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factored moment expansion in Model F. 

In summary, we have shown that the relativistic effects of nucleon recoil, Fermi motion 
and a reduced effective mass M* can be accounted for accurately in both Models F and 
G. Keeping only terms through O(p^) gives ~ 5% accuracy, which is small on the scale 
of the dominant kinematical effects. Keeping terms through O(p^) gives ~ 1% accuracy. 
We emphasize that in both Models F and G, the RCSR was evaluated using the same 
moment expansion coefficients (those in Appendix B), but in Model F the effective charges 
and magnetic moments depend on q and M* . This simplification requires using a factored 
moment expansion, in which Fermi motion effects entering though the effective form factors 
G^^, G\j^ are ignored. In the cases studied this is an excellent approximation, as explained 
further below. 



4.3 Discussion of Results 

One can understand the results just presented by looking individually at the electric and 



magnetic recoil functions te and tm-, defined in (|3.6| ) and (|3.7|) , as functions of M* . We 
shall look at the variations of e*(q) and /i*(q) of (3.14) in Model F, compared to e^- and 
in Model G. We are especially interested in the factored moment expansion in Model F. 
We shall examine the reasons for its success, and under what conditions such an approach 
may be expected to work for other off-shell assumptions. We argue that a factored moment 
expansion, like that adopted here, should be accurate in a large class of off-shell form factor 
models. 

In Figures 5a and 5b, we show the electric and magnetic recoil functions: r'£;(p, q) of 
( p.6| ) and (p, q) of ( p.7| ), for two values of M* . In both figures, the solid curve represents 
the electric term r£;(0,q), which dominates the Coulomb response for a system of Dirac 
nucleons. The magnetic term rA/(0, q) = 0, as can be seen from the second equality in ( |3.7| ). 
Thus the Coulomb response for a nucleon at rest is purely electric, as noted in Ref. [Q. 
The electric function r£;(0,q) decreases from unity at g = to its limiting value 1/2 as 
g — > oo. These limits are independent of M*, however, for smaller M*, r£;(0,q) decreases 
more rapidly with g, since it is a monotonic function of q/M* . This increased suppression is 
easily accommodated by using an appropriate value for M*, as discussed for the Dirac case 
in Figure 1. 

The corrections for Fermi motion to r^; and tm are also shown in Figures 5a and 5b, 
and are denoted by (5r£;(q) and 5r^/(q), respectively. These corrections were calculated by 
averaging r^(p, q) and rA/(p, q) over the Fermi sphere, taking pp = lA2 fm~^. For a system 
of Dirac nucleons, 6rE and 6rM enter the response with equal weighting. We noted in Ref. M 
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for M*/M= 1 that this leads to a nearly exact cancellation over the entire momentum range, 
and ensures that the lowest-order sum rule, RCSR-I, is accurate to within ~1%. Note that 
this cancellation persists to ~1% accuracy for M*/M<1, although both 6rE and Svm nearly 
triple for M7M~l/2, as in the MFT example. 

For nucleons with anomalous magnetic moments, this near cancellation of 5r£;(q) and 
(5rAf(q) no longer occurs, since these function are multiplied by different factors: e^, /i^ of 
(3.8) for Model G, or /if (q) of (3.14) for Model F. In Model G, (5rA/(q) is multiplied 

by the constant factor + ~ 11.4, while ep + e^ = l. For M*/M = l, this gives a 10-20% 
increase in the Coulomb sum S(q), as was discussed in Ref. 0. As M*/M decreases, the 
enhanced effect of 5rM(q) increases further: for M*/M = 0.556, the Fermi motion correction 
increases with q about as quickly as r£;(0,q) decreases, as shown in Fig. 5b. This leads to 
the behavior of S(q) for this value of M*, shown in Fig. 3b, and remarked on earlier, i.e., 
S(q) ~ 1. As M* is decreased to the MFT value, this increase in Stm effectively triples the 
importance of higher moments. This accounts for the decreased accuracy of the RCSR at 
C(p^), compared to that applied to free nucleons. 

Model F has qualitatively similar behavior, as we saw in Figs. 4a and 4b. Here, however, 
S(q) is modified by the effective charges and moments (3.14), which are shown in Fig. 6 
(solid curves). For this case {fip^+fJ'n'^) is somewhat smaller than (/ip+/i^) and increases with 
q over the range shown. This explains the smaller values of S(q) in Figs. 4a and 4b than in 
Figs. 3a and 3b, and the upturn with increasing q in Fig. 4b. 

We have seen that the convergence properties of the RCSR in Model F are similar to 
those in Model G, in spite of the use of a factored moment expansion that ignores certain 
Fermi motion effects, which would enter through the effective form factors (3.11). In order 
to understand this result, we first examine the effective form factors themselves. In Figure 6, 
we show these in the form in which they enter the recoil function ro-(p, q; Q"^) in Model F, i.e., 
[G*Ea/GEp? and [Glj^/GEpf- We have used (3.12) and M7M = 0.556 in their evaluation. 
The dashed curves show the form factors at p = 0, and the solid curves were obtained by 
averaging over the Fermi sphere with pp = 1.42 fm~^. (Note that the electric form factors 
tend to their free values in the limit >0, while the magnetic form factors tend to their free 
values in the limit g — oo. This can be seen from the definitions given in ( |3.1lD .) We see 
that the Fermi motion corrections to the ratios [G*/Gep]'^ are small compared to the ratios 
themselves. Thus terms O(p^) and higher are small, especially in the magnetic form factors. 
(Terms 0{p ■ q) vanish upon angle-integration with spherical symmetry, and therefore do 
not contribute to the results in Figure 6.) 

To understand the success of the factored moment expansion, consider the second-order 
sum rule, Sjj{q). It omits two types of Fermi corrections: 1) 0{p^) terms which come from 
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a moment expansion of the effective form factors and multiply r£;(0, q) and rM{0, q), and 2) 
products of linear terms involving (p ■ q), one term from the effective form factors, and one 
from the recoil factors ve and vm- The first of these make at most ~ 5% correction to the 
electric contribution, and less than ~ 1% correction to the magnetic, as can be seen from 
Figure 6. These terms can therefore be neglected to a good approximation. Linear terms do 
not enter the magnetic contribution, as can be seen from the (p x q)^ factor in ( p.TI) . Linear 
terms have been omitted from the electric contribution, and could in principle have been 
appreciable. In fact, these terms taken together result in less than a 1% error in the RCSR, 
as can be seen by comparing the dashed and dot-dashed curves for S'//(q) in Figures 4a and 
4b. Thus the success of the factored moment expansion in Model F can be attributed to the 
smallness of Fermi corrections in the effective form factors themselves, and in part to the 
complete absence of linear terms in the magnetic contribution. 

We believe that the success of the factored moment expansion in Model F may also 
occur for a much wider variety of off-shell models than we have considered here. The basis 
for this claim is that Fermi corrections tend to be small in general, representing effects 
of order (pi?/M*)^. The main exception is the magnetic correction, which is substantial 
because of two properties: rAf(0,q) = 0, and {G\j^ + G%j^) / G e-^ ~ 10. The magnetic 
contribution Stm is small, but is then enhanced by roughly an order of magnitude by the 
anomalous magnetic moments. No such enhancement occurs for 5rE- Thus in Model F, the 
further corrections due to p-dependence in [G*e^ / G e^i]'^ and [C^^^/G^ip]^ are negligible. In 
other off-shell models, for which ^'^(Pjq) is expressible in the form (3.9), we expect Fermi 
motion effects in the corresponding effective form factors also to be small. We therefore 
propose that a factored moment expansion, like that employed here in Model F, may also be 
accurate for the evaluation of Fermi motion effects a wide variety of off-shell models. This 
allows one to evaluate the RCSR using the expansion coefficients given in Appendix B. (We 
include 5r£;(q), although a small contribution in nucleon models with anomalous magnetic 
moments, to ensure the cancellation of Fermi corrections for the Dirac system, and to allow 
the "consistent" calculation of ^//(q) in Figs. 5a and 5b.) 

5 Summary and Conclusions 

The main results of this paper can be summarized as follows: The effects of nuclear binding 
in the mean-field approximation enter the Coulomb response function Wc{oj,c^ and the 
Coulomb sum S(q) in a way that can be characterized by a reduced effective mass M* . These 
effects depend on the behavior of the electromagnetic current for the off-shell kinematics 
in the nuclear medium. The sensitivity of the Coulomb sum S(q) to the choice of off- 
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shell model has been discussed previously |T0| , |Tl|]. We consider two illustrative models: 
F and G. We demonstrate in these models that the RCSR of Ref. can be extended to 
account for relativistic binding effects, in addition to the purely kinematic effects of recoil 
and Fermi motion treated in Ref. 0. The resulting RCSR is no longer model-independent, 
since one must make some assumptions about M* and the off-shell behavior of the form 
factors. However, to the extent that the dominant effects of the nuclear medium can be 
characterized by the simple parameter M*, the Coulomb sum rule analysis can be applied 
to data, with the goal of looking for nuclear structure effects beyond the mean field, e.g., 
two-body correlations. We emphasize that the form of the sum rule is sufficiently general 
to accommodate a broad class of off-shell form factor models, i.e., those for which the recoil 
function ro-(p, q) can be expressed in the form (3.9). 

The real change in transforming to the M* basis enters through the nucleons-only ap- 
proximation, under which the sum rule is derived. Strong potentials in relativistic models 
mix free N components into the initial and final interacting nuclear states, and modify the 
spacelike response function. Chinn, Picklesimer and Van Orden |]^, ^] have isolated the ef- 
fects of this mixing of N components, and show results qualitatively similar to those seen 
in Fig. 4. (Their form factors are closer to our Model F than G.) The transformation to 
an M* basis automatically incorporates these potential effects in a convenient representa- 
tion, although the use of an effective mass independent of position or momentum may only 
approximate the physical situation. 

How might one apply this RCSR to experimental data? First, the Coulomb sum S(q) 
must be calculated from the measured Coulomb response function, as in (|1.1|) . Then, one 
must make some specific assumptions about the off-shell behavior of the current operator, as 
discussed in Section 3. This should be cast in the G* form as in (3.10) and (3.11) for Model 
F (or G* = G{Q^) for Model G), with the substitution (3.12) to remove (approximately) any 
residual w-dependence from G%^/Gep and G\j^/Gep- The choice of model is not restricted to 
those presented in Sections 3.1 and 3.2. With the factored moment expansion, the relativistic 
recoil function (3.1) takes the form (3.13) and the ratios G* /Gep are evaluated at p = 0, as 
in (3.14). The recoil factor r^(q) is then evaluated in a moment expansion, e.g., to O(p^) 
as in (4.4), with (p^)cr fixed by other experimental information, or by a model as in (4.5), 
or as a free parameter. The expansion coefficients rjo^(q) are given in Appendix B, and are 
functions of the parameter M*. 

The modified Coulomb sum S'A(q) is obtained by forming the ratio (4.2) of the experimen- 
tally determined numerator S(q), and the (model dependent) recoil denominator ZrA(q). 
The expected behavior of S(q) with increasing q is that it will approach the one-body term 
S(^)(q) of (1.3), assuming that both C(q) and S[^^)(q)^0 in (1.2), as q^oo. The modified 
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sum S'yi(q) will then approach unity if the assumed form of the current is correct, and if the 
effective mass M* is appropriately chosen. Should that be the case for a given set of data, 
it would be reasonable to assume that the recoil functions have been correctly chosen. The 
sum rule (1.3) can then be used to investigate the two-nucleon correlation function in the 
ratio form 

5^(q)^l + C^(q) + -^, (5.1) 

by looking for deviations from unity at moderate momentum transfers. The two-body cor- 
relation function, 

Ca{cO ^ (5.2) 
Z r^(q) 

is related to the standard nonrelativistic correlation function. (See Eqs. (5.8) and (5.16) in 
Ref. P].) An interesting application of the RCSR is the constraint of the viable off-shell form 
factor models by analyzing the experimentally determined Coulomb sum S(q) in different 
off-shell models and comparing to S'A(q) — 1 beyond the expected range of correlations. 

A further remark about the analysis of data with the RCSR seems appropriate. It 
has become customary for the experimental Coulomb response data to be integrated in a 
modified form of (1.1), as suggested by DeForest |T^, in which the proton electric form factor 



Gep{Q^) is replaced by GEpiQ^) = GEpiQ^),J (l + r) / {1 + 2t), with r^Q^AM^ (see, e.g., 
Refs. |]14|, |T5[)- We explained in Ref. why this procedure will not lead to a non-energy- 
weighted sum rule: to obtain a sum rule of the form (1.2) or (1.5), the extra tu-dependence 
should not be introduced into the definition of the Coulomb sum. The kinematic effects 
included in the recoil factors r^(q) depend on the 3- momentum transfer q, rather than on 
the invariant Q^. It is this property which preserves the sum rule in passing from (1.2) to 
(5.1). 

The use of the M*-basis for both initial and final states, with the same value of M*, 
implicitly assumes that the nucleus is large enough that the nucleon kinematics in the final 
state are sensitive to the scalar field. For smaller nuclei, it may be necessary to account for 
the effects of finite size. This could possibly be accomplished through a rederivation of the 
RCSR in a Hartree representation based on bound, localized nuclear states. Another issue 



is the momentum-dependence of the mean-field potentials, as discussed in Refs. [|^, p^ , 
for example. Since the RCSR has been derived in momentum space, such a modification is 
relatively straightforward. As a first approximation, which is consistent with our conclusion 
that Fermi motion effects tend to be small, it seems reasonable to use M* (0) for initial states 
and M*(q) for final states. This prescription preserves the basic structure of the RCSR, in 
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that no further dependence on p is introduced. This would require the initial- and final-state 
masses to be treated distinctly, however, and would complicate the form of the coefficients 
in Appendix B. 

We are currently investigating the effects of virtual NN pairs (of mass M*) on the RCSR, 
as they enter intermediate excited states in the random phase approximation. These were 



considered previously by Horowitz |T^, and appear to be significant. We are also interested 
in how energy dependent terms, which can enter in off-shell models where a substitution of 
the form (3.12) is not appropriate, affect the RCSR. 
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Appendix A: in Terms of Sachs Form Factors 

An alternative form of the current operator (2.7) may be obtained by making a Gordon 
transformation EU| on matrix elements between free plane-wave spinors: 



1 



1 + r 



(Al) 



where P,i = l{p+p')ti, and 



(A2) 



and the F's and G's are related by (2.8). Another form of may be obtained |ST]: 



(A3) 



For /i = we have 



Po = -^{Ep+^ + Ep) , 



ro = 2757 ■ (p X q) 



(A4) 



For matrix elements between mean-field spinors, we have for /i = 0: 
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1 + f* 



2n 757 ■ (pxq) 



(A5) 



2M* ' '"'■^ ' 2M*^ 
In Model G, we use the on-shell forms of GEiQ"^) and GMiQ"^), given in (|A5|) . In Model 
F, we use G*^^ and G^/g., given in (3.11). Expression (A5) makes clear the origin of the 
functional forms in (3.6) and (3.7). 



Appendix B: Coefficients of Factored Moment Expan- 
sion 

In this appendix, we give the coefficients rjo-, for even powers through O(p^) in the 
moment expansion. These are written as they appear in Model G, i.e., in terms of the usual 
nucleon charge (e^) and magnetic moment (/Uo-): 



2 E*^ + M* 



2E* 



?^2a(q) 



q 

n*2 ;\ /f *3 



■51 



12E*^M*^ 





'E*^-M*' 


+ A 


3E*M*2 _ 



(Bl) 
(B2) 



240E*"M*^ E^ + M* 



GAE^" + 29E*^M* 



-QE*^M*^ + QE*^M*^ + l^E*^M*^ + 24E*^M*^ 



+30E*^M*^ - 90^*^M*^ - 210E*A'r 



-QAE*^ - 2AE*^M* + 16E*^M*^ 



(B3) 



-IQE*^M*^ - A^E*^M*^ - 2AEl^M*^ 



These forms are also applicable to Model F, upon making the replacements ta^e*^{c\) and 
/^a^K(q) using (3.14). 
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Figure Captions 



Fig. 1. Coulomb sum S(q) for a uniform system of Dirac protons, at two values of M*. Two 
different evaluations of the lowest-order RCSR are shown as indicated. 

Fig. 2. RCSR evaluated for a uniform system of nucleons with anomalous magnetic mo- 
ments, for M*/M — 1. The Coulomb sum E(q) in (4.1) and Sa{<1) in (4.2) are shown as 
indicated. 

Fig. 3. Same as Fig. 2, but for M*/M < 1 in Model G. Results are shown for (a) M*/M ^0.71S 
and (b) M*/M ^0.556. 

Fig. 4. Same as Fig. 3, but in Model F. Note that Sii{q) is calculated in two ways: the 
"factored" (dashed) and "consistent" (similar dot-dashed) moment expansions. 

Fig. 5. Electric and magnetic recoil functions, r£;(p, q) and rM(p, q), in a convenient sepa- 
ration. Results are shown for (a) M*/M^l and (b) M7M = 0.556. The functions rB(0, q) 
(solid), 5r£;(q) (dashed) and 5rM(q) (dot-dashed) are shown; rM(O,q)=0. (Note the scale 
changes for 5rE and 5rM-) 

Fig. 6. Effective form factors (squared) in Model F, for M*/M — 0.556. Solid curves are 
obtained by averaging (3.11) with (3.12) over the Fermi sphere, and include Fermi correc- 
tions to all orders. Dashed curves are e^^(q) and //^^(q) of (3.14), as used in the "factored" 
moment expansion. 
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